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Abstract 

We derive the kinetic luminosity function for flat spectrum radio jets, using the empirical and theoretical 
scaling relation between jet power and radio core luminosity. The normalization for this relation is derived 
from a sample of flat spectrum cores in galaxy clusters with jet-driven X-ray cavities. The total integrated jet 
power at z = is Wtot ~ 3 x 10"'" ergs s 

energy density deposited by jets as etot ~ ^ a lu ergsivipc . dulu wtot 
luminosity sources. Comparing etot to the local black hole mass density pbh gives an average jet production 
efficiency of ejct — ejot / Pbuc^ w 3%. Since black hole mass is accreted mainly during high luminosity states, 
ejot is likely much higher during low luminosity states. 
Subject headings: galaxies: jets — black hole physics — accretion 



^ Mpc By integrating Wtot over red-shift, we determine the total 
i 2 X 10^* ergs Mpc^'^. Both Wtot and etot are dominated by low 



\. INTRODUCTION 

The M — a relation between black hole m ass and the ve- 
locity dispersion of the host galaxy's bulge (iGebhardt et al.l 
|200C| Ferrarese & Merritt 2000) shows that the growth of 
black holes and large scale structure is intimately linked. X- 
ray observations of galaxy clusters show that black holes de- 
posit large amounts of energy into their environ ment in re- 
sponse to radiative losses of the cluster gas (e.g. Birz an et"aLl 
|2004). Finally, mechanical feedback from black holes is be- 
lieved to be respo nsible for halting star formation in massive 
elliptical galaxies (* Springel et al.ll2065l) . 

These arguments hinge on the unknown efficiency ejct with 
which growing black holes convert accreted rest mass into jet 
power. Constraints on ejct are vital for all models of black 
hole feedback. Efficiencies of 1 % are typically assumed, but 
this number is derived for phases of powerful jet outbursts 
only and neglects the trend for black holes to become more 
radio loud at lower accretion rates (Ho & Ulvestad'j 2001l: 
iGallo et al. 2003; Merloni et al. 2003; Falcke et al. 200|). 

The discovery of a tight correlation between core radio 
and X-ray lumino sity in accreting black hole X-ray binaries 
jGallo et al.ll2003h in low luminosity states inspired a number 
of theoretical investigations of how jet radio emission relates 
to the accretion state and ra te of the black hole. In the classi- 
cal model by IB landlord & Koenigll (1 19791) . the flat spectrum 
radio synchrotron emission of a compact jet core is produced 
by superposition of self-absorbed synchrotron spectra, each 
from a different region in the jet. The model predicts a de- 
pendence of the radio luminosity on jet power Wjct of 

the form oc W^JJ'^^. More generally. iHeinz & SunyaevI 
(l2003h showed that any scale invariant jet model producing a 
powerlaw synchrotron spectrum with index must obey the 
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Given a measurement of L^,, we can thus estimate a jet's ki- 
netic power, up to a multiplicative constant Wq (which we 
determine in m.2i . In f|2] we will use this relation to con- 
struct the kinetic jet luminosity function from the observed 
flat spectrum radio luminosity function ^l{L,^) (abbreviated 
as FSLF below). In Sj3] we will derive the current mean jet 
power per cubic Mpc and the jet production efficiency ejct of 
black holes. Section |4] summarizes our results. Throughout 
the paper we will use concordance cosmological parameters 
ofilu = 0.3, flA = 0.7, Ho = 75kms-iMpc-i. 

2. THE AGN-JET KINETIC LUMINOSITY FUNCTION 

As pointed out in iMerlonil (|2004 andlH einz et all (|2004 , 
we can use eq. ([T]) and the observed FSLF to derive the under- 
lying kinetic luminosity function $^ of flat spectrum jets: 



= $1 



12Wo\WoJ 

We will follow (iDunlop & Peacocklll99dl DP90) in using a 
broken powerlaw to describe ^^(Ly): 



$L(i.) = 



Lc{z) 



(3) 



c{z) J ' \Lc{z 
From DP90, we adopt ai — 1.85 and 02 ~ 3. oi is well de- 
termined at low z, but at higher z, the flux limit of the DP90 
sample approaches Lc and an accurate determination of ai 
is not possible anymore. In fact, within th e anti-hierarchical 
scenario for SMBH growth (e.g. Merloni 2004) a change in 
slope at high redshifts is expected, as more powerful black 
holes were more common at high redshift, as in deed observed 
in X- ray and optically selected AGN samples dHopkins et alj 
I2OO6I) . For lack of better information, we will assume a con- 
stant ai below. It is reasonable and, in fact, necessary that 
$L has a low-luminosity cutoff at some minimum luminosity 
imin (see fl221». 

The red-shift dependence of Lc and pq{z) varies with cos- 
mology. DP90 adopted a = 1, Ha = cosmology 
with Hq ~ 50kms^^ Mpc"'^, which gave -Z^o, (1,0,50) — 

6.7x1032 ergs Hz-is-ixl0^-^^[^"(^+^^"' '1. DP90 param- 
eterized po{z) as Po,(i,o,50) = 0.43 Mpc-^ x lol^- 
with cn = {-7.87, -5.74, 93.06, -738.9, 2248, -2399'}. 
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2.1. Correcting for relativistic boosting 

Since jets are relativistic, the observed FSLF is affected 
by Doppler boosting. The Dopppler-correction of luminos- 
ity functions has b een discussed i n a n u mber of publica- 
tions, most notably lUrrv & ShafeJ (11984 : lUrrv & Padovanil 
( 1199 lb . Following these authors, we will neglect the contri- 
bution from the receding jet. The error introduced by this 
approximation is small compared to the other sources of un- 
certainty. The Doppler factor for a jet with Lorentz factor 
r = — p^), velocity (3 = v/c and viewing angle ^los 

is then given hy S = l/[r (1 — /?cos (0los))]> with a maxi- 
mum of (5i„ax = 1/(1 + — The jet luminosity for a 
flat spectrum source is then boosted by a factor 6^. 

Without knowledge of the underlying jet four- velocity dis- 
tribution, exact Doppler correction is impossible. However, 
for sensible distributions that show a clear peak at some 
Tmcan, it is Sufficient to approximate the distribution as a delta 
function, allowing decomposition. We will assume that the 
velocity distribution is well behaved in such a way. 

A rest-frame (i.e., intrinsic) FSLF of the form of eq. (|3]l, 
subject to Doppler boosting, will still be observed as a bro- 
ken powerlaw with the same indices, but with an additiona l 
powerlaw regime with slope of -3/2 (lUrrv & ShafedfT984h . 
At low z, the low luminosity slope of the observed FSLF is 
well determined to be steeper than —3/2, indicating that the 
—3/2 Doppler tail must lie at luminosities below the DP90 
flux limit. There are indications of a turnover to a —3/2 slope 
at luminosities below Lmin.obs ~ 10^^ eiKsHz"^ in the 



low red-shift, low luminosity sample of ("Nagar et al.l l2005i 
NFW), indicating that the rest frame FSLF becomes shallower 
than -3/2 below Lmin = imin.obs/^max- ^igh z, fli is not 
well enough determined to dr a w this conclusion. 

Foflowing lUrrv & ShafeJ (11984 . the rest frame (i.e., 
Doppler corrected) FSLF must be of the form of eq. (O, with 
the observed break luminosity ic,obs and the low-luminosity 
cutoff/turnover £min.obs each Doppler boosted by a factor 



and the normalization corrected by a factor 



A 



[(1 - /3)3-2''i - 1] 

/3(2ai - 3) 



(4) 



2.2. Estimating the kinetic power normalization Wq 



To estimate the n ormalization Wp o f the radio — :j et-power 
relation in eq. ([T]l, ' Heinz et aTl ( 12004 and iHeinz & Grimml 
(|2005) used information from three well studied radio galax- 
ies, M87, Cygnus A, and Perseus A, for which estimates of 
the kinetic power from large scales (and kpc jet-scales in the 
case of M87) exist, along with measured flat spectrum fluxes 
for the jet core. 

Given the recent X-ra y surveys of ga l axy clusters with c en- 
ti-al radio sources (Birza n et al] 120041; lAllen et all |2006), a 
more sophisticated estimate of Wq is now possible. Taking 
all radio sources with robust kinetic power estimates based 
on X-ray cavities and with observed nuclear flat spectrum 
fluxes, we composed a sample of 15 sources. Whenever more 
than one measurement of the jet kinetic power was available 
we have taken the logarithmic average of the available data. 
Fig. [Uplots the core flux against the kinetic power for the 13 
sources considered. Other estimators for jet power are avail- 
able in the literature (e.g. Willott et al. 1998), based on total 
steep spectrum radio power. Since, unlike cavity based mea- 
surements of W, these estimates are rather model dependent, 
we will not employ them here. 
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Fig. 1 . — Time averaged kinetic power plotted against flat spectrum core 
radio power. Lines show the best fit normalization (with one sigma uncer- 
tainties). 

Eq. (|T) predicts that the two quantities Wq and should 
be related by a power-law with index 12/17. To derive the 
constant of proportionality VFo we performed a least-squares 
fit to the data in Fig.[T] fixing the slope at 12/17. The best fit 
for Wq is shown as a thick solid line, along with two lines for 
the 1 -sigma uncertainty derived from the scatter in the plot. 
Abitrarily fixing Lq at 



Lo,obs = 7 X lO^^ergsHz 



(5) 



and introducing the parameter wj^^ = Vl^o/10*"'ergss ^ for 
convenience, the best fit for Wq is 



= itli X 10'''^ ergs s"^ = waa x 10'^'' ergs 



(6) 



The observed core luminosity Lo.obs must be corrected for 
Doppler boosting. For an unbiased sample of sources, the 
core luminosity will, on average, be de-boosted by a factor 
(5|q = 4/[r2(2 — /3)^], corresponding to an average view- 
ing angle of 60°. However, we cannot assess what selection 
biases affect the sample of sources contributing to this esti- 
mate (e.g.. X-ray sensitivity to detect jet-indu ced cavities se- 
lects flga/nif beamed sources : ,EnBlin & Heinzt2002 ). We will 
write Lo — io,obs/^60 witii the implicit understanding that 
(Jgo subsumes the unknown effects of any line-of-sight bias in 
the sample from Fig.[T] 

After taking all corrections into account, the un-boosted ki- 
netic luminosity function is 



17po.o 
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where we defined the critical power Wc as 



T A2 

-'^c.obs^gO 

C'O.obs'^njax 



(7) 



(8) 



has a low luminosity slope of (^{W) oc H^(5-i7ai)/i2^ 
This implies that the total integrated kinetic power is dom- 
inated by the lowest power sources for values of ai steeper 
than ai > 29/17 1.7. Given that DP90 measured ai w 
1.85, the luminosity function must have a break or cutoff at 
some somewhere below the flux limit of DP90 (possi- 
bly given by the value of imin,obs indicated by NFW), since 
otherwise the total kinetic luminosity would diverge. 



The jet production efficiency of growing black holes 
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Fig. 2. — Integrated jet power density for different parameter combinations 
at z = 0, as a function of mean jet Lorentz factor Fmean- Solid lines: 
ai = 1.85, dotted lines: ai = 1.8, dashed lines: ai = 1.7, dash-dotted 



2.5 X lO^'^ergsHz- 
1 (NFW). 



s-i (DP90), 



lines: ai = 1.6; black lines: Lmin = 
grey lines: Lmin = lO^'^ ergs Hz~^ i 

3. DISCUSSION 

The total power released by jets from all sources under the 
FSLF (i.e., all black holes except those contributing to the 
observed steep-spectrum luminosity function, abbreviated as 
SSLF), per comoving Mpc^^, is simply the first moment of the 
kinetic luminosity function: 



Wtotiz) 



dx 



llpQ.ohs{z) Wc 



(9) 



where Xo{z) = (£min.obs(z)/ic,obs(2)) For Lmin ~ 

10^'^ ergsHz^i s^^ (NFW), this yields 

Wtotiz = 0)«2.8 X 10*°ergss-iMpc-3 2?it;44 (10) 

24 

^60 \ 1 
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2-2ai 
max 



(11) 



This is the estimated total power released per Mpc^ by flat 
spectrum jets today. The main sources of uncertainty in Wtot 
are Wo, ai, Fnioan, and, for large values of ai, the estimate 
of Lrnin- Fig. 12] shows Wtot/it'44 as a function of Fmcan for 
different values of imin and ai. 

To put this value of Wtot in context, it is useful to com- 
pare it to the average stellar luminosity density L^, 2 x 
10^^ ergs Mpc^'^ (Ellis et al. 1996) and to th e current su- 



perno va power VFsN ~ 10'^''' ergs s^^ Mpc~^ dMadau et alj 
119981) . which is well below the integrated jet power. 

The estimate of Wtot is dominated by low luminosity 
sources (which is why imin is critical for Wtot)- This sug- 
gests that AGN make a much bigger contribution to feedback 
in regular galaxies than commonly assumed. It also implies 
that low luminosity AGN dominate the global kinetic energy 
output of black holes in a quasi-steady state, rather than short, 
energetic bursts of individual black holes. This picture agrees 
well with the concept of slow, "effervescent" feedback envi- 
sioned to be responsible for AGN heat ing in galaxy clusters 
(lBegelmanll200UlChurazov et al.ll2002l) . 

In order to derive the total energy density released by jets, 
we can integrate Wtot from eq. (|9]l over redshift, taking the 
proper cosmological corrections into account. We will use 
an upper limit of Zmax — 5 for the redshift integral, but the 
results are not sensitive to the exact value of ^max- Taking the 
most conservative approach by using the observed radio flux 



limit from DP90 as a solid upper limit on the low-luminosity 
cutoff Lmin, WC find a lower limit of 



etot 



dt. 



dz—Wtot > 1.5 X 10'"ergsMpc"-'Pw44 
dz 



(12) 

on the total integrated jet energy density. A more realistic 
assumption would be to adopt the low-redshift value of imin 
from NFW and assume the same redshift evolution for imin 
as that of L^. This gives an estimated value of 

etot ~ 1-7 X lO^^ergs Mpc"^ I? u;44 



(13) 



Keeping Lmin fixed at the z = value from NFW at all 
redshifts provides a robust upper limit of etot < 2.5 x 

j^gSS ergs M-pC-'-^VWii. 

Comparing etot to the mean cosmi c black hole mass dens ity 
of PBH « 3.3 X 10^ Mq Mpc-3 (lYu & Tremainel2002h fi- 
nally yields the average conversion efficiency of accreted rest 
mass to jet power for supermassive black holes: 



ejct 



etot 
Pbhc^ 



3%X'W44 



(14) 



for the same assuptions that went into eq. ( fT3] l. Assuming 
that Ljnin{z) is smaller than the survey flux limit and larger 
or equal to Lmin.NFW- gives limits of 0.25% 2? ^44 < ejot < 
4.3% 15^44. This estimate of ejot is broadly consistent with 
previous best-guess estimates of the jet conversion efficiency, 
typically believed to be of the order of 1% - 10%. 

Note, however, that the estimate of pen from 
lYu & Tremaind (l2002l) includes mass accreted in all phases of 
black hole growth. Black holes grow predominantly through 
radiatively efficient accretion: The observed amount of X-ray 
background radiation is equivalent to about 10% of the total 
black hole rest mass energy measured today, implying that, 
for typical radiative efficiencies of order 10%, most of the 
accreted mass must have contributed to the product ion of 
the X -ray background in a radiatively efficient mode ( SoltanI 
|1982|) . Radiatively efficient accretion flows are typically 
radio quiet (i.e., inefficient at producing jets). Thus, the 
conversion efficiency during low-luminosity accretion phases 
must be significantly higher than the average e implied by 
eq. (fT4l ). Given that about 10% of AGN are radio loud, the 
average black hole accumulates a fraction of /,„ — 90%/9o 
of its mass during radio quiet, radiatively efficient accretion. 
Thus, the jet conversion efficiency ejot during radio loud 
phases must be at least a factor of (1 — fm)^^ ^ 10 larger 
than shown in eq. ( fT4b . of order 30%. 

Efficiencies of several tens of percent are not implausible, 
however: If black hole sp in is important in jet launching 
dBlandford & Znaie3ll977 ') and if black holes accrete large 
amounts of angular momentum (i.e., are close to maximally 
rotating for a significant f raction of their life, a s suggested by 
recent merger-tree models IVolonteri et al.ll2005i) . they can lib- 
erate up to 30% of the accreted rest mass energy by black hole 
spin extraction alone, implying ejot of up to 42%. 

Finally, we will briefly discuss the possible contribution 
of steep-spectrum sources to ejct- By definition, the FSLF 
contains all black holes except those included in the SSLF 
(which are dominated by optically thin synchrotron emission). 
The scaling relation from eq. ([T]l does not hold for steep- 
spectrum sources. We can, however, derive an upper limit 
on the contribution from steep spectrum sources. Assuming 
a typical optically thin synchrotron spectral index of 0.65, 
any underlying flat spectrum component would have to fall 
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below 20% of the observed 5GHz luminosity, otherwise the 
source would become too flat to qualify as a steep-spectrum 
source. Using the SSLF from DP90, the same redshift integral 
that yielded eq. ( fT2b provides an upper limit on the contribu- 
tion from the flat spectrum sources underneath the SSLF of 
Cstccp < 6% 2? W44, where T) is the Doppler correction for the 
steeper spectral index. Since the SSLF is shallow at low lu- 
minosities, Estccp is dominated by sources around Lc and not 
affect by the uncertainty in imin- 

4. SUMMARY 

Starting from the relation between kinetic jet power and flat 
spectrum core radio luminosity, we derived the kinetic lumi- 



nosity function of flat spectrum radio sources. We found that 
the kinetic luminosity density is dominated by the lowest lu- 
minosity sources, indicating that constant, low level efferves- 
cent type heating is important in black hole feedback. Inte- 
grating the kinetic luminosity density over redshift and com- 
paring it to the estimate of the current black hole mass density 
showed that the efficiency of jet production by black holes is 
of the order of a few percent and smaller than 10%. However, 
since most of the power comes from low luminosity sources, 
which are not believed to contriubte much to the total mass 
accretion of the black hole, the efficiency of jet production 
during low luminosity, jet-driven phases must be significantly 
higher. 
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